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Abstract. We report measurements of the real and imaginary parts of the AC conduc- 
tivity in the quantum limit, hu > ksT of insulating nominally uncompensated n-type 
silicon. The observed frequency dependence shows evidence for a crossover from inter- 
acting Coulomb glass-like behavior at lower energies to non-interacting Fermi glass-like 
behavior at higher energies across a broad doping range. The crossover is sharper than 
predicted and cannot be described by any existing theories. Despite this, the measured 
crossover energy can be compared to the theoretically predicted Coulomb interaction 
energy and reasonable estimates of the localization length obtained from it. Based on a 
comparison with the amorphous semiconductor NbSi, we obtain a general classification 
scheme for electrodynamics of electron glasses. 



1. Introduction 

Strong electronic interactions are known to play a central role in disordered 
solids, of which Coulomb glasses are a canonical example. The lack of 
metallic screening on the insulating side of the metal-insulator transition 
(MIT) enables long-range Coulomb interactions [1]. Efros and Shklovskii 
(ES), following the original considerations for the non-interacting Fermi 
glass case of Mott[2], derived a form for the T = K photon assisted 
frequency dependent conductivity describing the crossover from interacting 
Coulomb glass- like behavior to Fermi glass-like behavior [3]. These deriva- 
tions were based on a theory of resonant absorption[4] and take into account 
the mean Coulomb interaction between two sites forming a resonant pair 
U(r u ) = e 2 /sir UJ , where = t;{ln(2Io/hLu)] is the most probable hop 
distance between pairs and £\ is the dielectric constant. The real part of 
the ES crossover form for the frequency dependent conductivity is: 

o-i = /3e 2 g 2 fio[ln(2I /hio)} 4 [hLO + U(r w )\ (1) 

where j3 is a constant of order one, go is the non-interacting single 
particle density of states (DOS), Iq is the pre-factor of the overlap integral 
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and £ is the localization length. The concentration dependent localization 



where x is the dopant concentration, x c is the critical dopant concentration 
of the MIT (x c = 3.52 x 10 18 cm~ 3 in Si:P[5]) and v is the localization length 
exponent. 

Neglecting logarithmic factors, Eq. (1) predicts a gradual crossover from 
linear to quadratic behavior as the incident photon energy exceeds the 
interaction energy of a typical charge excitation. For the case where the 
photon energy, hw > U{r w ), one recovers the quadratic frequency depen- 
dence, plus logarithmic corrections, that Mott originally derived for the 
non-interacting Fermi glass case[2]. In the opposite limit, fouj < U(r u ) the 
conductivity shows an approximately linear dependence on frequency, plus 
logarithmic corrections, and the material is called a Coulomb glass. We 
should note that Eq. (1) was derived for the case where Hu > A, the 
Coulomb gap width. However a quasi-linear dependence (albeit with a 
different pre- factor) and an eventual crossover to Mott's non-interacting 
quadratic law is still expected even for the case where ftw < A. 



There is a lack of experimental evidence to either corroborate or dis- 
prove Eq. (1) due to the difficulties associated with performing frequency 
dependent measurements in the so-called quantum limit, i.e. Hlo > ksT, but 
at small enough photon energies so as to not be exciting charge carriers to 
the conduction band. Moreover, in order to study the possible crossover 
from Mott to ES type behavior, one must measure across a broad enough 
bandwidth centered about the characteristic crossover energy scale for in- 
stance the Coulomb interaction energy U or the Coulomb gap width [6], 
A. 

There have been some very recent experiments that have attempted to 
address these issues. M. Lee et al. found that for concentrations close to the 
MIT the expected linear to quadratic crossover occurs, but is much sharper 
than predicted [7] . They proposed that this sharp crossover was controlled 
not by the average interaction strength U as in Eq.(l) [3], but instead by 
a sharp feature in the density of states, i.e. the Coulomb gap [8]. They 
postulated that this Coulomb gap was not the single particle one measured 
in tunneling, but rather a smaller "dressed" or renormalized Coulomb gap 
that governs transport. There is some evidence from DC transport that 
such a feature exists, at least close to the MIT [9]. 




approached 



1 UJE\ 



(2) 



Wln(2I /hu) 
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2. Experiment 

Nominally uncompensated n-type silicon samples were obtained from Recti- 
con Enterprises Inc. A Czochralski method grown boule with a phosphorous 
gradient along its length was cut into 1 mm thick discs. Room temperature 
resistivity was measured using an ADE 6035 gauge and the dopant concen- 
tration calibrated using the Thurber scale[10]. The Si:P samples discussed 
here span a range from 39% to 69%, stated as a percentage of the sample's 
dopant concentration to the critical concentration at the MIT. A number of 
samples were measured before and after etching with a 4%H F +96%i? NO3 
solution; this resulted in no difference in the results. 

In the millimeter spectral range, 80 GHz to 1000 GHz, backward wave 
oscillators (BWO) were employed as coherent sources in a transmission 
configurationfll]. The transmitted power through the Si:P samples as a 
function of frequency was recorded. For plane waves normally incident 
on a material, resonances occur whenever the thickness of the material 
is an integer number of half wavelengths. Both components of the com- 
plex conductivity can be uniquely determined for each resonance. The real 
part of the conductivity was evaluated at microwave frequencies from the 
measured loss of highly sensitive resonant cavities at 35 and 60 GHz via 
the perturbation method. This is a common technique and is described in 
the literature[12]. The conductivity as determined from the resonant cavity 
data was normalized to the DC conductivity at higher temperatures, at 
above approximately 25 K. The resonant cavity data confirmed the linear 
dependence on frequency of the real part of the complex conductivity into 
the microwave regime for the samples closest to critical. 

3. Results 

In Fig. 1, we show the T^0 frequency dependent conductivity for two 
samples. This data, representative of all samples in our range, shows an ap- 
proximately linear dependence at low frequencies and then a sharp crossover 
to an approximately quadratic behavior at higher frequencies. This is the 
qualitatively expected behavior from Eq. (1). However, as seen by the over- 
layed fits, Eq. (1) provides only a rough guide. The solid lines are linear and 
quadratic fits to the low frequency and high frequency data respectively. 
The dotted line is a fit to the form of the ES crossover function achieved 
by summing the separately determined linear and quadratic fits. As can 
be seen, the crossover between linear and quadratic portions is much more 
abrupt than the ES function predicts. The dashed line is a fit using the same 
method as Ref. [7], namely forcing the linear portion to pass through the 
low frequency data, as well as the origin and leaving the pre-factor of the 
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Figure 1. Low temperature frequency dependent conductivity data plotted as a function 
of frequency. The Si:P samples shown are at 50% and 69% dopant concentration relative 
to the critical concentration. The solid lines are linear and quadratic fits to the lower and 
upper portions of the data respectively. The dotted and dashed lines are fits following 
the form of Eq. (1) using two different methods described in the text. 

quadratic term as a free variable. The fit is not satisfactory in either case. 
A sharp crossover as such is observed over our entire doping range and has 
been observed previously in an analogous system, Si:B, for samples closer 
to the MIT [7]. Note, that a linear dependence is seen in the imaginary 
part of the conductivity o<i over the whole measured frequency and doping 
range [13]. This is consistent with theoretical predictions [14]. 

Because our data spans a large range of concentrations, the doping 
dependence of the crossover energy scale can be analyzed to see whether 
its dependence is consistent with other energy scales, e.g. the Coulomb 
interaction energy U or the Coulomb gap width A as per Ref. [7]. Recall 

that the Coulomb interaction energy between two sites forming a resonant 

2 

pair is U{r^) = which is dependent on concentration via the dielectric 
constant (measured, but not shown) and the localization length dependent 
most probable hop distance. By equating the crossover energy scale to the 
expected functional form for this Coulomb interaction energy we are able to 
determine the magnitude of the localization length and its scaling exponent. 
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With an appropriate pre-factor in the overlap integral [15], Io = 10 13 s _1 
for the expression for the most probable hop distance term, r u , we get a 
localization length dependence of £ oc (1 — x/x c ) -0 ' 83 with a magnitude 
of 21.2, 19.9, 20.1, 14.5, 14.3 and 13.0 nm for the 69%, 62%, 56%, 50%, 
45% and 39% samples respectively. The localization length exponent is 
close to unity, the value originally predicted by McMillan in his scaling 
theory of the MIT [16], and the magnitude of the localization length is 
reasonable. Due to the fact that we obtain reasonable estimates for the 
relevant physical parameters over the whole doping range, we do not favor 
the previous speculation that it is in fact the Coulomb gap energy that 
creates the sharp crossover and hence sets its energy scale [7]. 

The approximately linear power law seen in the Coulomb glass regime 
at low lo in Fig. 1 of the conductivity can be expressed with the imagi- 
nary part as a simple Kramer-Kronig compatible form, (t{oj) = A{iu) a = 
Ato a cos + iAuj a sin (^p) ■ In order to determine the power a one can 
take the ratio of |cr 2 1 versus o\ (with the frequency as a variable). The 
power a is given by, 



Fig. 2a shows the ratio mentioned above of the imaginary to the real 
part of the dielectric constant for Si:P. Similar data from amorphous NbSi 
is included for comparison purposesfl]. We note that this ratio for Si:P 
is large and essentially constant across the entire doping range. From Eq. 

(2) , one expects \a\\ to be approximately equal to |cr 2 1 to within a factor 
of 2-5 (with a reasonable estimate for Io) because |cr 2 1 oc E\ ■ uj. Applied 
to Si:P, the theory correctly predict a linear correspondence between o\ 
and | (72 1) but incorrectly predicts the proportionality by at least a factor 
of thirty. The proportionality is near the predicted value for NbSi, but 
has a dependence on the doping concentration which is presumably due 
to entering the quantum critical (QC) regime as discussed below. Here 
we have used the susceptibility 4irx of the dopant electrons (i.e. with the 
background dielectric constant of silicon subtracted) in the expression for 
the magnitude of the imaginary component of the conductivity in Eq. 3. 

The middle panel in Fig. 2 shows the power a as determined by Eq. 

(3) . The values for Si:P are approximately equal to, but slightly less than 
one, consistent with Fig. 1. This indicates that the prefactor of the real and 
imaginary components of the complex conductivity have the same concen- 
tration dependence. The situation is different for NbSi. Near the MIT, a{u) 
is expected to cross over to the QC dynamics [17, 18], i.e. o\ oc uj 1 / 2 when 
£, the localization length, is of the same scale as £ u , the dephasing length 
(the characteristic frequency dependent length scale) [19]. This should be 




(3) 



ArmitageCG.tex; 2/02/2008; 1:20; p. 5 



6 



10' 



« 10 



3 



1.0 - 



0.8 



10' 



* * 1 ! 



• Si:P 
□ NbSi 



0.3 0.4 0.5 0.6 0.7 
1 - X/X 



Figure 2. The upper panel shows the ratios of the imaginary to the real part of the 
complex conductivity for samples of Si:P and amorphous NbSi. The NbSi data is adapted 
from Ref 1. The middle panel shows the calculated powers of a as determined from Eq. 
3. The dashed line through the NbSi data is a guide to the eye. The bottom panel shows 
the divergence of the prefactor of the real part of the conductivity, and the dotted line 
is a simple power law fit. 



a smooth crossover and therefore looking at a fixed window of frequencies, 
a continuous change from ui — > lo 1 ^ 2 is expected, similar to that measured 
for NbSi shown in the middle panel of Fig 2. 

Setting the relations for localization length and dephasing length equal 
[19], one finds the crossover condition for the frequency in terms of the 

normalized concentration, uj cx £q ^1 — j-^j where z is the dynamic expo- 
nent. As the prefactor can vary from system to system, the fact that we see 
an a ss 1 across our entire doping range in Si:P, but an a that approaches 
1/2 in NbSi indicates that the critical regime in Si:P is much narrower and 
out of our experimental window. This is consistent with simple dimensional 
arguments [20] that show the crossover should be inversely proportional to 
the dopant density of states. The much smaller dopant density in Si:P vs. 
NbSi (a factor of 10 3 ) is consistent with a narrower QC regime as compared 
to NbSi. 

The bottom panel in Fig. 2 shows magnitude of the real part of the 
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conductivity as the MIT is approached. This demonstrates that the pref- 
actor A can be written as a function of the normalized concentration, i.e. 
A oc /(l - ^) for Si:P. 

4. Discussion 

In typical interacting systems, the effects of correlations become simpler 
as one goes to lower energies and/or lower temperatures. The canonical 
example of this is a Fermi liquid where at T=0 and co=0 one recovers the 
non-interacting theory, but with parameters that are substantially mod- 
ified (renormalized) from the free electron ones. The Coulomb glass is a 
fundamentally important example in solid state physics because it belongs 
to a class of systems where this does not occur and the non-interacting 
functional forms are not recovered at asymptotically low energies. 

As predicted in the crossover function Eq. 1, we have observed in Si:P 
that the non-interacting functional form is recovered in the high-frequency 
limit and the low-frequency response shows interactions. Within the theory 
this is a result of the additional internal excitation structure of a resonant 
pair caused by interactions and the fact that any one pair can be thought 
of as a distinct entity, well separated in energy from other spatially nearby 
pairs because of the large disorder induced energy spread. This internal 
structure enables the excitation of pairs relatively deep within the Fermi 
sea and changes the factor in the initial state phase space from uj to U (r) + 
uj. In contrast, in amorphous NbSi a smooth crossover is observed as the 
MIT is approached from a linear frequency dependence to a power law 
characterized by an exponent smaller than unity. This is consistent with the 
eventual frequency dependence a oc uj 1 / 2 expected from quantum critical 
scaling arguments. No crossover to uj 2 was observed in NbSi. 

Since it is the exponent a that distills the important physics (it indicates 
the phase space of initial states for a = 1 or 2 and the presence of critical 
dynamics for a = 1/2), we propose that one can classify the electrody- 
namics of electron glasses based on their a value. A schematic showing the 
parameter space for a is shown in Fig. 3. Here one has a's close to 1/2 
near the MIT. There is a smooth crossover to Coulomb glass- like a = 1 at 
lower doping levels and an intervening non-interacting Fermi glass regime 
at even higher energies and lower dopings. We expect that these general 
considerations are valid, despite the fact that some of the parameter space 
may not be accessible in certain systems. For instance, one may begin to 
excite structural modes at energies high enough to see uj 2 behavior in NbSi. 
For dopings not close enough to critical in Si:P, excitation to the conduction 
band may be observed before critical dynamics are. 
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Figure 3. A schematic showing the parameter space for values taken by a for a oc ui a . 
Note that the boundaries drawn on the plot are smooth crossovers and not sharp onsets. 
A classification based on a gives a taxonomy for the electrodynamics of electron glasses. 

5. Conclusion 

In summary, we have observed a crossover in the frequency dependence 
of the conductivity from Coulomb glass-like behavior to Fermi glass-like 
behavior across our entire range of doping concentrations in Si:P. The 
existence of a crossover is consistent with theoretical predictions, but it 
is sharper than predicted. The fact that we see the same functional form 
over the whole doping range (even deep into the insulating regime, where 
Eq. 1 is expected to be more valid) shows that the nature of the low energy 
charge excitations is qualitatively the same over the whole doping range; 
the inadequacy of Eq. 1 in describing the frequency dependent conduc- 
tivity quantitatively is not limited to concentrations close to critical. In 
the amorphous semiconductor NbSi, we observe a gradual crossover from 
a = 1 Coulomb glass-like behavior to quantum critical-like dynamics. This 
comparison allows us to obtain a general classification scheme for the elec- 
trodynamics of electron glasses based on the exponent of the frequency 
dependence a. We expect that this classification or 'taxonomy' will be valid 
even when certain regimes are not experimentally accessible. 
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